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I 

O 1 ABSTRACT 

^2 ' Context. The distance to the Galactic Centre (GC) is of importance for the distance scale in the Universe. The value derived by Eisenhauer et 
al. (2005) of 7.62 ± 0.32 kpc based on the orbit of one star around the central black hole is shorter than most other distance estimates based on 
a variety of different methods. 

Aims. To establish an independent distance to the GC with high accuracy. To this end Population-II Cepheids are used that have been 
^ ' discovered in the OGLE-II and OGLE-III surveys. 

Methods. Thirty-nine Population-II Cepheids have been monitored with the SOFI infrared camera on 4 nights spanning 14 days, obtaining 
typically between 5 and 1 1 epochs of data. Light curves have been fitted using the known periods from the OGLE data to determine the mean 
i\-band magnitude with an accuracy of 0.01-0.02 mag. It so happens that 37 RR Lyrae stars are in the field-of-view of the observations and 
mean if -band magnitudes are derived for this sample as well. 
Results. After correction for reddening, the period-luminosity relation of Population-II Cepheids in the if-band is determined, and the derived 
slope of —2.24 ± 0.14 is consistent with the value derived by Matsunaga et al. (2006). Fixing the slope to their more accurate value results 
in a zero point, and implies a distance modulus to the GC of 14.51 ± 0.12, with an additional systematic uncertainty of 0.07 mag. Similarly, 
from the RR Lyrae if -band period-luminosity relation we derive a value of 14.48 ± 0.17 (random) ± 0.07 (syst.). The two independent 
determinations are averaged to find 14.50 ± 0.10 (random) ± 0.07 (syst.), or 7.94 ± 0.37 ± 0.26 kpc. The absolute magnitude scale of the 
adopted period-luminosity relations is tied to an LMC distance modulus of 18.50 ± 0.07. 
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1. Introduction 

The distance to astronomical objects is a crucial parameter, yet 
often very difficult to obtain with high precision. The distance 
to the Galactic Centre (GC) is of special importance, e.g for 
dynamics (Oort constants, determining distances using a rota- 
tion model), or for calibrating standard candles. The classically 
accepted value comes from the review by Reid (1993) and is 
R = 8.0 ± 0.5 kpc. 

Over the last few years the distance to the GC based on the 
orbit of the star called S2 around the central black-hole (BH) 
has caught attention. Initially, Eisenhauer et al. (2003) derived 
a value of 7.94 ± 0.42 kpc which was revised by Eisenhauer 
et al. (2005) to 7.62 ± 0.32 kpc having more epochs of data 
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available. The neglect of post-Newtonian physics in these anal- 
ysis may have lead to an underestimate of the distance by about 
0.11 ± 0.02 kpc (Zucker et al. 2006), leading to a current best 
estimate of 7.73 ± 0.32 kpc (corresponding to a distance mod- 
ulus (DM) of 14.44 ± 0.09) to the GC based on the BH. 

On the other hand, most other recent distance determi- 
nations give a longer distance, more in line with the classi- 
cal value: (1) High-amplitude delta-scuti stars give 7.9 ± 0.3 
kpc (McNamara et al. 2000); (2) RR Lyrae stars suggest a 
value of 8.8 ± 0.3 kpc (Collinge et al. 2006), 8.3 ± 1.0 kpc 
(Carney et al. 1995) or 8.0 ± 0.65 kpc (Fernley et al. 1987). 
(3) Earlier work on the Red Clump gave a longer distance of 
8.4 ± 0.4 kpc (Paczyriski & Stanek 1998), although Nishiyama 
et al. (2006) derive 7.52 ± 0.10 (stat) ± 0.35 (syst) kpc, and 
Babusiaux & Gilmore (2005) 7.7 ± 0.15 kpc; (4) From a com- 
parison of Miras found in the OGLE database in the direc- 
tion of the Galactic Bulge (GB) to those in the Magellanic 
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Clouds, Groenewegen & Blommaert (2005) find a distance in 
the range 8.5 to 9.0 kpc, in agreement with earlier work on 
Miras (Catchpole et al. 1999); (5) Analysis of the Hipparcos 
proper motions of 220 Cepheids lead to i?o = 8.5 ± 0.5 kpc 
(Feast & Whitelock 1997); (6) Modelling the observed colour- 
magnitude diagram in V, I and J, K using a population synthe- 
sis code, Vanhollebeke et al. (2008) derive a distance of 8.60 ± 
0.16 kpc. 

With the exception of some Red Clump based distances, 
the results obtained by Eisenhauer et al. imply a much shorter 
distance to the GC than found by most other methods, and this 
calls for an independent investigation of this matter. 

In this paper the distance to the GC is determined us- 
ing Population II Cepheids (hereafter P2C) discovered in the 
OGLE micro-lensing survey, and for which the mean K- 
band magnitude will be determined by infrared monitoring. 
Comparing to the calibrated P2C period-luminosity (PL) re- 
lation in the /f-band from Matsunaga et al. (2006; hereafter 
M06) then provides the distance, after correction for redden- 
ing. 

In addition, the mean A'-band magnitude will be deter- 
mined for RR Lyrae stars that are in the field, and compared to 
the calibrated X-band PL-relation from Sollima et al. (2006). 
The Matsunaga et al. and Sollima et al. relations both imply an 
LMC DM of 18.50 as detailed in Sect. 5. 

In Sect. 2 the sample is discussed, and the observations are 
presented in Sect. 3 for the P2C and Sect. 3 for the RR Lyrae. 
The results are discussed in Sect. 5. 

2. The sample of Population n Cepheids 

Population-Il Cepheids are old, low-mass stars. They are the 
progeny of hot HB stars that after the exhaustion of core He- 
burning, move toward lower effective temperatures (post-early 
AGB), thus crossing the Cepheid instability strip. They are sys- 
tematically brighter than RR Lyrae stars and have periods rang- 
ing from slightly below one to a few tens of days 

Kubiak & Udalski (2003; hereafter KU) have searched the 
OGLE-II database for P2C Cepheids and found 54 objects. KU 
determined a period-luminosity relation in the reddening-free 
Wesenheit index (based on V, I photometry) and compared it 
to P2C Cepheids in the LMC. The difference in DM at a typ- 
ical period of logP = 0.5 is 3.58. By assuming a DM for the 
LMC of 18.5 (50.1 kpc) the distance to the GC P2C becomes 
9.6 kpc (with a substantial error bar of 1 .5 kpc). This places the 
objects nominally in the GB region. Additional evidence is that 
KU found these objects to be located in a bar (see their Fig. 2), 
like the RR Lyrae (Collinge et al. 2006), a strong indication 
that the P2C Cepheids indeed are physically part of the GC re- 
gion. These properties gave us the confidence that determining 
individual distances to these objects would indeed provide an 
accurate distance determination to the GC. 

A subset of stars was selected based on the following crite- 
ria: The /-band light curves of these stars were inspected and 
only "smooth" ones were retained that show no or little effect 
of shocks in the atmosphere (Figure [T] shows two examples). 
The range in Galactic longitude was largely restricted to -2.5 

^+2.5 degrees to avoid any additional smearing in distance 
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Fig. 1. Phased OGLE /-band light curves of two of the P2C. 
The periods of the two stars are 0.76522 (bul_sc24_1609) and 
9.94431 (bul_sc39_5385) days. 



due to the effect of the inclination of the bar. The entire period 
range from 0.8 to ±10 days should be covered. OGLE /-band 
images were inspected for crowding. 

In a later stage, it was also possible to extent the sample 
with the first results from the ongoing OGLE-III survey and 
a sample of 70 new P2C in the longitude range -2.5 to +2.5 
degrees could be added. After inspecting the light curves and 
/-band images of these stars as well, a final combined sample 
of 49 stars remained, of which 39 were actually monitored. 

The basic properties of the sample are listed in Table Q] in 
the first columns: Identification, R.A., Declination, pulsation 
period, galactic coordinates. 

3. The observations 

The observations were carried out with the SOFI infra-red cam- 
era on the 3.5m NTT on ESO/La Silla in the nights of 2007, 
June 24, 28, July 3, 8 in visitor mode. 

Photometric conditions were excellent on the second and 
third night with seeing as low as 0.6" and the telescope was 
actually defocused. Weather conditions were poorer on the first 
and last night with seeing in the range 1.5-2" and cirrus and 
thin clouds. As the measurements of the P2C will be relative 
to the 2MASS objects in the field this additional extinction did 
not influence this program. 

Typically it was tried to observe the longest period (>7 d) 
Cepheids at the beginning and end of the night, the intermediate 
period ones (>3 d) three times per night, and the shorter period 
ones were observed one after the other over the entire night. In 
total we obtained 362 epochs of data of 39 P2C. 
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Images in the K s band (hereafter simply Jf-band) were 
taken with the shortest possible on-chip integration time of 
DIT= 1.2 seconds using a pixel-scale of 0.288" and resulting 
in a field-of-view of almost 5x5 arcmirQ. The "auto-jitter" ob- 
serving block (OB) was used with 9-13 exposures. The relative 
large number was used to have sufficient redundancy in creat- 
ing a sky image in these relatively crowded fields. The if-band 
was chosen for the observations as the dispersion in the PL- 
relation is smaller in that band than in J or H (M06). 

The data were reduced with the newly released SOFI data 
reduction pipelin^. The pipeline takes into account the cross- 
talk, the flat-field (dome flats taken with a special OB), bad 
pixel cleaning, and images correlation and reconstruction. 

The reduced images were trimmed to the original 1024 x 
1024 pixel size to eliminate the under exposed edges. The as- 
trometric solution was done using the WCSTools suit43 match- 
ing the stars in the field against the 2MASS catalog. Typically 
several tens of 2MASS objects are available and the rms in the 
solution typically is 0.2". Source extraction and PSF photom- 
etry was done using version 1.3 of DoPhot (Schechter et al. 
1993). A dedicated Fortran program was written to match the 
2MASS objects with the sources in the field, and determine 
the offset between the instrumental magnitude and the 2MASS 
if s -band magnitude. The brightest and faintest 2MASS objects 
were excluded to avoid problems of saturation in the SOFI im- 
ages and poor quality 2MASS data. A straight line was fitted 
to the data, and this determined the offset to be applied to the 
SOFI instrumental magnitudes to put them on the 2MASS sys- 
tem (an example is shown in Fig. |2j. The typical error in the 
determination of this offset is 0.006 mag, and this is added in 
quadrature to the error in the instrumental magnitude for each 
star to give the total error in the observed magnitude. 

This error is verified by comparing the magnitudes over the 
different epochs for field stars and estimating the precision with 
which the mean magnitude can be determined. 

After having determined the photometric offsets, the light 
curves of the P2C could be constructed. The light curves were 
then fitted with a sine-function with the known period from the 
OGLE data, and this allowed the determination of the mean 
X-magnitude with high precision. Figure[3]shows the (phased) 
light curve for two P2C Cepheids with a short and long period 
in the sample. The last three columns of Tab. Q] list the mean 
X-band magnitude, the precision in this determination and the 
number of epochs. 

4. RR Lyrae stars 

Collinge et al. (2006) present a catalog of 1888 fundamental- 
mode RR Lyrae stars in the OGLE fields. It was verified which 
of them happen to be located in the SOFI fields of view, and 
for those the if -band light curve was extracted, and fitted with 
a sine-curve with the known period. Table|2]present the results, 
with columns 2-4 taken from Collinge et al. Column 5 lists 
in which P2C field it is located and in the analysis the corre- 
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Fig. 2. The difference between SOFI instrumental magnitude 
and 2MASS magnitude versus instrumental magnitude for a 
single epoch for two different objects. The symbols indicate 
different 1 magnitude bins in 2MASS magnitude. 



sponding reddening for that field is taken (see below). The last 
columns lists the mean magnitude, the error, and the number of 
epochs. Figure|4]shows the light curve for 2 stars. 

5. Analysis 

5.1. P2C: K -band PL-relation 

The if -band magnitudes are de-reddened using the model by 
Marshall et al. (2006). They present 3-dimensional i\T s -band 
extinction along 64 000 lines of sight at a resolution of 15' in 
typically four distance bins. Using VizieR0 the available data 
within 20' radius of the targets were retrieved. The mean and 
dispersion was determined of the Ak values in the bin that cor- 
respond to a distance larger than 4 kpc (for 2 stars a smaller 
distance had to be adopted). Table Q] lists in Columns 7-9 the 
mean, dispersion and the number of data points used. 

Columns 10 of this table list the Ak value derived by mul- 
tiplying by 0.12 (based on the reddening law of Cardelli et al. 
1989) the Ay value for the OGLE-II GB fields by Sumi (2004). 

Based on the MACHO survey towards the GB Popowski et 
al. (2003) derive visual extinction for 9717 elements at a reso- 
lution of about 4'. The available data within 5' radius of the tar- 
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Fig. 3. The X-band light curves of two P2C. The top panel shows the magnitude versus time, the bottom panel the phased light 
curve. 
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Fig. 4. The If -band light curves of two RR Lyrae. The top panel shows the magnitude versus time, the bottom panel the phased 
light curve. 
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Name 




RA 


DEC 




(J2000) 


(J2000) 


bul24.1609 


17 


53 44.41 


-32 57 14.0 


bul2 14.2.135 


17 


57 35.12 


-28 26 11.6 


bul214.6_53717 


17 


56 15.76 


-28 14 02.3 


bul39.2239 


17 


56 05.04 


-29 54 51.8 


bul22.3993 


17 


56 33.03 


-30 36 33.5 


bul3 1.662 


18 


01 56.92 


-28 55 11.6 


bul43.351 


17 


35 08.13 


-27 31 25.8 


bull63.7_46320 


17 


51 54.32 


-31 56 39.8 


bul30.1107 


18 


01 47.33 


-29 07 39.1 


bul215.7_107864 


17 


58 54.55 


-28 21 30.8 


bul3_791 


17 


53 27.82 


-30 19 55.3 


bul5_3719 


17 


50 21.49 


-29 54 27.5 


bul32.2167 


18 


02 56.05 


-28 40 39.2 


bul4_170 


17 


54 38.21 


-30 10 41.8 


bul44.5324 


17 


49 30.95 


-29 50 58.2 


bull88. 1.11087 


17 


59 30.71 


-30 19 59.7 


bul2 1.3035 


17 


59 53.52 


-28 56 37.5 


bul45.1189 


18 


03 33.61 


-30 01 14.6 


bul23.611 


17 


58 14.58 


-31 33 23.7 


bul34.4631 


17 


58 01.53 


-28 59 56.6 


bul20.961 


17 


58 55.93 


-29 10 57.0 


bul41.3841 


17 


51 42.51 


-32 41 41.3 


bul4_8846 


17 


54 06.61 


-29 16 22.1 


bul4_2323 


17 


54 55.52 


-29 57 31.0 


bul39.616 


17 


55 12.35 


-30 07 24.1 


bull 95. 4.144166 


17 


53 40.70 


-28 56 47.0 


bul22.815 


17 


56 46.47 


-31 07 07.3 


bul3_1755 


17 


53 34.75 


-30 12 39.7 


bul38.3260 


18 


01 32.68 


-29 49 11.5 


bull80.5_151030 


17 


51 57.55 


-30 21 39.2 


bul20.3867 


17 


58 55.95 


-28 43 19.5 


bul2_1657 


18 


04 19.87 


-29 02 44.2 


bul25.693 


17 


54 24.38 


-33 06 51.0 


bul40.2524 


17 


50 58.92 


-33 08 53.1 


bul214.5_169880 


17 


56 37.75 


-28 01 28.4 


bul7_540 


18 


08 44.36 


-32 13 10.6 


bul39.5385 


17 


55 23.12 


-29 31 35.3 


bul334.7_59737 


17 


40 50.79 


-25 03 52.7 


bul333.2_25311 


17 


36 34.05 


-27 18 08.3 



Period 
(days) 



I 



b Ai 



A 



A 



W 



K 



0.76522 
0.82684 
0.83088 
0.92213 
0.93110 
0.93951 
1.09776 
1.21123 
1.33916 
1.40015 
1.48412 
1.50105 
1.50523 
1.53185 
1.55140 
1.69825 
1.73020 
1.74795 
1.85470 
2.33010 
2.88445 
3.20082 
3.51336 
3.54254 
3.65017 
3.82791 
4.49944 
4.57875 
5.53097 
6.89131 
7.13700 
7.45600 
7.73934 
7.77000 
8.07957 
9.52109 
9.94431 
11.38492 
14.89102 



357.44 
1.77 
1.80 
0.33 

359.77 
1.82 

359.97 

358.12 
1.63 
1.98 

359.68 

359.70 
2.14 

359.94 

359.66 
0.33 
1.58 
1.04 

359.13 
1.33 
1.27 

357.45 
0.66 
0.16 
0.05 
0.90 

359.36 

359.80 
1.00 

359.49 
1.67 
1.97 

357.37 

356.98 
2.02 

359.64 
0.58 
2.74 
0.33 



-3.56 
-2.00 
-1.65 
-2.46 
-2.89 
-3.07 
2.71 
-2.71 
-3.14 
-2.22 
-2.18 
-1.39 
-3.14 
-2.32 
-1.20 
-3.31 
-2.69 
-3.92 
-3.68 
-2.37 
-2.63 
-3.06 
-1.76 
-2.26 
-2.40 
-1.52 
-3.19 
-2.14 
-3.44 
-1.91 
-2.40 
-3.59 
-3.76 
-3.16 
-1.61 
-5.95 
-2.13 
2.94 
2.56 



0.284 
0.317 
0.352 
0.275 
0.330 
0.193 
0.492 
0.433 
0.215 
0.282 
0.323 
0.310 
0.178 
0.278 
0.313 
0.271 
0.198 
0.177 
0.312 
0.234 
0.223 
0.346 
0.295 
0.271 
0.268 
0.315 
0.270 
0.318 
0.214 
0.381 
0.233 
0.174 
0.275 
0.344 
0.367 
0.149 
0.278 
0.456 
0.449 



0.018 
0.028 
0.005 
0.026 
0.051 
0.019 
0.045 
0.078 
0.016 
0.035 
0.049 
0.023 
0.012 
0.027 
0.039 
0.045 
0.019 
0.019 
0.009 
0.031 
0.038 
0.028 
0.025 
0.022 
0.023 
0.024 
0.031 
0.046 
0.027 
0.042 
0.031 
0.007 
0.020 
0.022 
0.050 
0.004 
0.025 
0.036 
0.034 



0.302 



0.316 
0.329 
0.217 
0.440 

0.229 

0.347 
0.688 
0.193 
0.311 
0.720 



0.220 
0.324 
0.302 
0.233 
0.318 
0.311 
0.311 
0.316 

0.329 
0.347 
0.220 

0.233 
0.186 
0.281 
0.353 

0.160 
0.316 



0.354 0.011 6 
0.343 0.023 6 
0.237 0.021 3 



0.233 0.015 4 
0.278 0.011 4 



0.212 0.011 5 



0.219 0.014 5 

0.169 0.009 5 

0.312 0.033 5 

0.219 0.004 5 



0.286 0.010 5 



0.206 0.018 4 

0.243 0.012 5 
0.189 0.014 4 



0.462 0.017 5 
0.323 0.006 4 



14.237 
13.805 
14.118 
13.836 
13.583 
13.925 
13.629 
12.573 
12.126 
13.088 
14.255 
13.304 
13.351 
13.027 
13.449 
12.904 
13.159 
13.125 
12.783 
12.690 
12.661 
12.196 
12.449 
13.267 
12.281 
12.101 
12.415 
12.613 
12.175 
11.601 
11.639 
11.464 
11.775 
11.999 
12.180 
11.220 
11.086 
14.200 
10.856 



0.010 
0.049 
0.057 
0.023 
0.013 
0.010 
0.007 
0.010 
0.006 
0.009 
0.022 
0.028 
0.008 
0.010 
0.011 
0.018 
0.017 
0.018 
0.007 
0.024 
0.006 
0.008 
0.012 
0.112 
0.007 
0.014 
0.035 
0.055 
0.121 
0.007 
0.007 
0.007 
0.006 
0.006 
0.007 
0.006 
0.007 
0.018 
0.008 



13 
5 



13 
11 
11 
11 

8 
7 
11 
10 

8 
9 
9 
9 
10 
9 
9 
9 



Notes: 

(1) Reddening from Marshall et al. (2006). (2) Reddening based on Sumi (2004). (3) Reddening based on Popowski et al. (2003). 



gets was retrieved. The mean and dispersion was determined, 
and multiplied by 0.12 to obtain the if -band extinction, which 
is listed with the number of elements used as Columns 11-13 
in Table □ 

The comparison of the reddenings based on the Marshall et 
al. model and the values from Sumi and Popowski et al. suggest 
that there are no systematic effects and that the error in the 
adopted reddening is realistic. 

The observed A"-magnitudes are then de-reddened, and 
Fig. [5] shows for the P2C the observed period-luminosity rela- 
tion in the A' -band, together with the best-fitting line (exclud- 
ing the cross), K = (-2.24 ±0.14) (log P - 1.2) + (10.578 ± 
0.099), with an rms of 0.41 mag. This is based on a weighted 
least-squares fit, where the error in Ak is added in quadrature 



to the error in K. Eliminating the possible three outliers near a 
period of 1.3 days gives a slope of —2.41 ± 0.08 and a ZP of 
10.529 ± 0.059 with an rms of 0.28 mag. 

This dispersion is largely due to the intrinsic depth of 
the Bulge which is of the order of 1 kpc (e.g. Babusiaux & 
Gilmore 2005), and was already seen in the dispersion around 
the Wesenheit Pi-relation by KU. 

Plotting the residuals versus galactic longitude and latitude 
suggested a slight dependence on these parameters, and a gen- 
eral linear fit was made (excluding only one object): 

K = (-2.24 ± 0.13)(logP - 1.2) + (10.60 ± 0.17) 

+ (-0.028 ± 0.031) I + (0.005 ± 0.031) b (1) 
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Table 2. The sample of RR Lyrae stars, ordered by period. 



Name 


RA 


DEC 


Period 


P2C field 


K 


a 


n remark 




(J2000) 


(J2000) 


(days) 











bul40.259253 

bul39.15393 

bul3 1.680405 

bul45.245148 

bul20.43721 

bul24.369931 

bul45.256663 

bul3_239578 

bul22.1 33009 

bul44.201914 

bul20.44361 

bul4_207615 

bul3 1.23420 

bul22.2 19909 

bul39.156007 

bul34.159387 

bul40.107313 

bul4_451052 

bul5 .261524 

bul22.402472 

bul4_439650 

bul32.87588 

bul22.144794 

bul25.187102 

bul22. 144296 

bul43.130470 

bul22.333396 

bul24.527005 

bul3_250926 

bul34.159754 

bul45.245339 

bul34.145186 

bul30.604974 

bul34.389919 

bul39.15358 

bul24.5 17291 

bul4_207388 



17.85087 
17.92076 
18.04767 
18.05646 
17.98309 
17.89280 
18.05978 
17.89180 
17.93951 
17.82289 
17.98202 
17.90976 
18.03285 
17.94488 
17.92021 
17.96703 
17.84697 
17.91352 
17.84023 
17.94676 
17.91289 
18.05029 
17.94068 
17.90468 
17.94163 
17.58358 
17.94312 
17.89388 
17.89106 
17.96594 
18.05951 
17.96658 
18.03144 
17.96842 
17.92118 
17.89567 
17.90973 



-33.1519 
-30.1590 
-28.6609 
-30.0464 
-29.1750 
-32.9691 
-30.0235 
-30.3038 
-30.6361 
-29.8616 
-29.1714 
-30.1619 
-28.9437 
-31.1046 
-29.5568 
-28.9742 
-33.1318 
-29.9600 
-29.8885 
-31.0980 
-29.9957 
-28.6886 
-30.6044 
-33.1508 
-30.5937 
-27.5280 
-30.5768 
-32.9262 
-30.2063 
-29.0003 
-30.0290 
-29.0164 
-29.1517 
-29.0303 
-30.1017 
-32.9539 
-30.1789 



0.39951 
0.43787 
0.43940 
0.44022 
0.45910 
0.47306 
0.47416 
0.48308 
0.48700 
0.48806 
0.48950 
0.50215 
0.50610 
0.52411 
0.52499 
0.54158 
0.54215 
0.54390 
0.55190 
0.55482 
0.55603 
0.57740 
0.57910 
0.57918 
0.58562 
0.59590 
0.61422 
0.61678 
0.63179 
0.63539 
0.63597 
0.66322 
0.67131 
0.67304 
0.69769 
0.76522 
0.78642 



bul40_2524 
bul39.616 
bul32_2167 
bul45_1189 
bul20.961 
bul24_1609 
bul45_1189 
bul3_791 
bul22_3993 
bul44_5324 
bul20.961 
bul4_170 
bul3 1.662 
bul22.815 
bul39.5385 
bul34_4631 
bul40_2524 
bul4_2323 
bul5_3719 
bul22.815 
bul4_2323 
bul32_2167 
bul22_3993 
bul25.693 
bul22_3993 
bul43.351 
bul22_3993 
bul24_1609 
bul3_1755 
bul34_4631 
bul45_1189 
bul34_4631 
bul30_1107 
bul34_4631 
bul39.616 
bul24_1609 
bul4_170 



14.899 
12.685 
14.536 
14.410 
14.399 
13.951 
14.421 
14.605 
14.331 
11.237 
14.994 
14.405 
12.923 
12.718 
14.119 
14.437 
14.430 
14.185 
14.320 
14.408 
14.419 
13.994 
14.364 
14.244 
14.261 
14.246 
14.525 
14.356 
13.017 
14.066 
14.034 
14.193 
14.361 
14.334 
14.323 
14.237 
14.090 



0.021 
0.017 
0.011 
0.010 
0.013 
0.008 
0.018 
0.082 
0.011 
0.006 
0.021 
0.028 
0.020 
0.007 
0.026 
0.017 
0.024 
0.013 
0.027 
0.014 
0.011 
0.014 
0.016 
0.015 
0.012 
0.013 
0.036 
0.011 
0.013 
0.016 
0.045 
0.025 
0.025 
0.016 
0.040 
0.010 
0.015 



5 
6 
11 



12 

7 
4 
9 
6 
6 
7 
5 
6 
4 
7 
5 



6 
7 
9 
9 
5 
9 

10 
4 

10 
6 
9 
6 
4 
8 
7 
5 

13 
6 



outlier 



outlier 



outlier 
outlier 



outlier 



The slope agrees within the error bar with that of the PL- 
relation derived by M06 which is also on the 2MASS system: 
M Ko = (-2.41 ± 0.05)(logP - 1.2) + (-4.00 ± 0.02) with 
a dispersion of 0. 14 mag. 

In order to derive the distance to the GC a Monte-Carlo 
simulation was performed. Data sets with new values for 
Ak and K were generated based on the mean values and 
(Gaussian) errors listed in Table [U and unweighted fits were 
made for a fixed slope of —2.41. This results in a zero point 
(ZP) of 10.512 ±0.013. 

The DM to the GC therefore is 14.5 1 ± 0.02, corresponding 
to 7.99 ± 0.09 kpc. The error takes into account the precision 
in our ZP and the error in the ZP of the M06 relation. 

Even adding in quadrature an additional 0.015 mag to the 
error on the photometry (allowing for an underestimate of the 
adapted error on the transformation from instrumental magni- 
tude to the 2MASS system), and doubling the error on the Ak 
value (with a minimum error of 0.03) increases the formal error 





1 ' 






_ □ 
m 
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0.5 


1 





Log (period) 

Fig. 5. The if-band PL-relation for P2C in the Bulge. The line 
is a best-fit excluding the cross. 



bar on the derived ZP only to 0.025 mag, and the error on the 
DM to 0.03 (0.12 kpc). 
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5.2. P2C: the Wesenheit-index PL-relation 

An alternative method to derive the distance is to use a 
reddening-free Wesenheit index. The most appropriate one is 
WIK = I — a(I — K) as the /-band magnitudes are available 
from OGLE. Unfortunately the available /-band data (Pritzl et 
al. 2003) for stars with if -band data in M06 is insufficient to 
derive an empirical relation, but one can use theoretical PL- 
relations as derived in several bands by Di Criscienzo (2007). 
These theoretical relations are in good agreement with ob- 
servations: The slope of the theoretical /-band PL-relation, 
— 2.10±0.06, is in agreement with the observed one of — 2.03± 
0.03 (Pritzl et al. 2003), while In the LT-band the theoretical 
slope of —2.38 ± 0.02 is in good agreement with the observed 
—2.41 ± 0.02 in M06. For an assumed mixing length parame- 
ter of 1.5, one needs to adopt a fiducial [Fe/H] value of —1.63 
in the theoretical PL-relation to obtain the observed / magni- 
tude at a typical period of log P = 0.5 (Pritzl et al. 2003), and 
similarly a fiducial [Fe/H] of —2.55 to obtain the observed K 
magnitude at log P = 0.5 (M06). 

The theoretical relation is, depending on the coefficient 
a in the Wesenheit-relation (Di Criscienzo, private comm.) 
WIK = (-1.13 ± 0.02) - (2.53 ± 0.04) logP + (0.06 ± 
0.01)[Fe/H] for ct\ = 1.252 (standard reddening law) or, 
(-1.17±0.02) - (2.57±0.04) logP+ (0.06±0.01)[Fe/H] for 
o?2 = 1-316 (allowing for anomalous reddening in /, see Sumi 
2004). For an [Fe/H] value of -2.0 ± 0.4 one expects ZPs of 
-1.25 ± 0.04 ( ai ) and -1.29 ± 0.04 (a 2 ) 

A potential complication is that the theoretical relation is 
derived for the Bessell & Brett (1988) photometric system, 
while in particular the OGLE /-filter is slightly non-standard. 
An appropriate transformation is applied to both the 2MASS 
K, and OGLE / to the Bessell & Brett system. It turns out 
however that for the typical (V — I) colours of the P2C in 
the observed fields and the particular value of the slope in the 
Wesenheit-index the corrections largely cancel, and the effect 
on the derived ZP is minimal. 

For a>i the derived slope from the observations is —2.25 ± 
0.17, for ct2 — 2.27±0.18. The agreement between the observed 
and theoretical slopes is less good (about 1.6-1.7(7 deviation) 
than that in the LT-band (only 1.3(7 deviation). 

When the slope is fixed to the theoretical value one ob- 
tains ZPs of 13.19 ± 0.01 and 13.04 ± 0.01, respectively, based 
on a Monte-Carlo simulation. Depending on the reddening law 
adopted one obtains a DM of about 14.44 (7.73 kpc), or 14.33 
(7.35 kpc) with an internal error bar of 0.04 mag (0.14 kpc). 

5.3. RR Lyrae 

Sollima et al. (2006) derive an empirical PLK-relation based on 
15 GC and the LMC cluster Reticulum, which for a metallicity 
of [Fe/H]= —1.0 (the average metallicity of RR Lyrae in the 
GB, see Walker & Terndrup 1991), reads M K = (-2.38 ± 
0.04) logP — (1.13 ± 0.13), which was calibrated against the 
trigonometric parallax of RR Lyra (Benedict et al. 2002). 

Removing five bright outliers the PL-relation becomes (cf. 
Figure [6}: 

Kq = (-1.36 ± 0.49) log P + (13.63 ± 0.17) 



lo : 




-0.4 -0.3 -0.2 -0.1 

Log (period) 

Fig. 6. The LT-band PL-relation for the RR Lyrae stars in the 
Bulge. The line is a best-fit excluding the outliers labeled in 
Table|2]which fall outside the plot. The dotted line is the fit for 
a fixed slope of —2.38. 

+ (+0.019 ± 0.021) I - (0.025 ± 0.043) b (2) 

The error in the derived slope is large but formally agrees 
within 2ct with the empirical slope by Sollima et al. For a fixed 
slope of —2.38, the ZP becomes 13.39 ± 0.13 (based on a 
Monte-Carlo simulation), resulting in a DM of 14.52 ±0.18 
(8.0 ± 0.7 kpc). 

6. Discussion 

The PL-relation in the LT-band of P2C in the GB is derived. 
The slope is found to be in agreement with that derived by M06. 
Fixing the slope to their more accurate value implies a DM to 
the GC of 14.51 with a formal error bar of 0.03. 

There is also the systematic error bar to consider. M06 pre- 
sented JHK period-luminosity relations based on 46 P2C with 
periods between 1.2 and 80 days in 26 Galactic globular clus- 
ters (GCs). For the absolute magnitude scale they adopted a re- 
lation between absolute V-magnitude of the Horizontal Branch 
and metallicity (Gratton et al. 2003), which in turn is calibrated 
using main-sequence fitting to three GCs. This calibration im- 
plies an RR Lyrae based LMC DM of 18.50 ± 0.09 (Gratton 
et al. 2003). M06 show that the DM based on P2C in the LMC 
and their A'-band PL-relation is also compatible with 18.5. 

They also show that there is no significant trend with 
metallicity over the range —2.2 Ss [Fe/H] Ss —0.5, in agree- 
ment with theoretical predictions (Bono et al. 1997, 2003, Di 
Criscienzo et al. 2007), and indicating that this LT-band PL- 
relation should be applicable for GC P2C as well. The metallic- 
ity of the P2C in the Bulge is unknown but that of RR Lyrae is 
estimated to be on average [Fe/H]= —1.0 (Walker & Terndrup 
1991). Any difference between that metallicity and the mean 
metallicity of about [Fe/H]= -1.5 of the GCs in M06 would 
result in an uncertainty in the ZP of <0.03 mag. 

Di Criscienzo et al. (2007) show that adopting a different 
My- [Fe/H] relation has a negligible effect on the derived slope 
of the NIR PL-relations. 

The ZP in the calibrating relation by Gratton et al. has a 
formal error of 0.07 and this has to be considered as a source 
of systematic uncertainty in the derived distance. 
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There is other source of (random) error to consider, namely 
how representative this particular set of 39 stars (minus the one 
outlier) that defines the PL-relation is in view of the fact that 
they scatter along the line-of-sight due to the intrinsic depth 
of the Bulge. To simulate this, additional Monte-Carlo sim- 
ulations were carried out. Random samples of 38 stars were 
selected from the original sample, and the Pi-relation re- 
derived. The dispersion in the ZP is about 0.11 mag. This is 
likely a slight overestimate as in this approach the randomly 
drawn samples do not necessarily have the large spread in pe- 
riod that the true sample was selected to have. 

The DM to the GC we derive based on the P2C is 14.51 ± 
0.12 (random) ± 0.07 (syst). The random error could be im- 
proved further by observing additional systems when the full 
OGLE-III database becomes available. 

Based on the serendipitously observed RR Lyrae stars in the 
field a DM of 14.52 ± 0.18 is derived. The error bar is for 50% 
due to the uncertainty in the adopted absolute magnitude of RR 
Lyra itself (Sollima et al.). Their PL-relation led to an LMC 
distance of 1 8 .54 ± 0. 1 5 . If instead we would assume the LMC 
distance to be 18.50 (to be consistent with the P2C calibration) 
then we would find a DM of 14.48 ±0.13 (random) ± 0.07 
(syst), were the systematic error comes from the uncertainty 
in the ZP of the observed LMC PL-relation. Like for the P2C 
sample, there is an additional 0.11 mag systematic uncertainty 
due to the limited sample size. The final DM to the GC based 
on the if -band RR Lyrae stars is 14.48 ± 0. 17 (random) ± 0.07 
(syst). 

As the two distance estimates are derived independently, 
they can be averaged and the best empirical estimate of the DM 
to the GC based on the current data is 14.50 ± 0.10 (random) 
± 0.07 (syst). 

The theoretical WIK relation gives a formal result of 14.44 
(or 14.33 with anomalous reddening) with an internal error bar 
of 0.04 mag. A random error of 0.11 has to be added to this, 
due to the limited sample, and as the theoretical relation is tied 
to the observed relations of M06 a similar systematic error bar 
of 0.07 has to be considered. 

Although within the error bar of the purely empirical re- 
sults, it brings up the question of reddening and the reddening 
law. An additional absorption in if of 0.1 mag would bring 
all three methods in very good agreement. On the other hand, 
the reddening estimates listed in Table 1 are in excellent agree- 
ment and are based on (J - K) colours (Marschall et al. 2006), 
(V - I) (Sumi 2004), and (V - R) (Popowski et al. 2003). If 
the absorption in K were underestimated, it would also imply 
an underestimate of the reddening in the other maps, or a sig- 
nificantly higher selective reddening A^/Ay ~ 0.16 instead 
of0.12. 

A final remark is that independent distances to some of 
these P2C may be obtained using surface-brightness relations 
(e.g. Groenewegen 2004) and the Baade-Wesselink technique. 
This would require better sampled if -band light curves than 
were needed for the present study and well-sampled radial 
velocity curves. Although observationally expensive it would 
give an improved understanding on the systematic error in the 
present analysis. 
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